There and (slowly) back again: 
Entropy-driven hysteresis in a model of DNA overstretching 
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When pulled along its axis, double-stranded DNA elongates abruptly at a force of about 65 
pN. Two physical pictures have been developed to describe this overstretched state. The first 
proposes that strong forces induce a phase transition to a molten state consisting of unhybridized 
single strands. The second picture instead introduces an elongated hybridized phase, called S- 
DNA, structurally and thermodynamically distinct from standard B-DNA. Little thermodynamic 
evidence exists to discriminate directly between these competing pictures. Here we show that within 
a microscopic model of DNA we can distinguish between the dynamics associated with each. In 
experiment, considerable hysteresis in a cycle of stretching and shortening develops as temperature is 
increased. Since there are few possible causes of hysteresis in a system whose extent is appreciable in 
only one dimension, such behavior offers a discriminating test of the two pictures of overstretching. 
Most experiments are performed upon nicked DNA, permitting the detachment ('unpeeling') of 
strands. We show that the long-wavelength progression of the unpeeled front generates hysteresis, 
the character of which agrees with experiment only if we assume the existence of S-DNA. We also 
show that internal melting (distinct from unpeeling) can generate hysteresis, the degree of which is 
strongly dependent upon the nonextensive loop entropy of single-stranded DNA. 

PACS numbers: 
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I. INTRODUCTION 

DNA in vivo experiences protein-mediated tensile 
forces large enough to alter the structure and stability of 
the hybridized state. Single-molecule manipulation tech- 
niques permit the study of these deformations through 
the controlled application of forces and torques [lj. In- 
triguingly, DNA pulled along its axis elongates suddenly, 
by a factor of about 1.7, at a force of 65 pN. Experi- 
mental thermodynamics does not conclusively establish 
the nature of this overstretched state. The abruptness 
of elongation suggests that the extended form represents 
a thermodynamic phase different from B-DNA. The pic- 
ture of overstretching as a force-induced 'melting' to un- 
hybrizided single strands @, H H H H is bolstered by 
the dependence of the overstretching force on parameters 
like pH and ionic strength [f| 0, [I[ ; these parameters are 
known to change the melting temperature of the molecule 
at zero tension. In addition, the binding of ethidium to 
DNA increases the overstretching force in a manner con- 
sistent with data from thermal melting studies [Ej]. The 
competing picture, which considers overstretching to be 
a transition to S-form DNA, a putative hybridized state 
1.7 times longer than B-DNA @, H ©, M, El El, is 
motivated by the mechanical properties of overstretched 
DNA. The latter has a stretch modulus of about 1600 
pN/0.34 nm Q, which exceeds the stretch moduli of both 
B-DNA (~ 1300 pN/0.34 nm) and unhybrizided forms of 
DNA at comparable forces. However, neither picture can 
be conclusively ruled out on the basis of these data. 

Here we show that these thermodynamic data can be 
used to develop a model of DNA able to discriminate be- 
tween the dynamics associated with the two pictures. We 



have constructed such a model, with a resolution of one 
basepair, whose kinetics we explore using a Monte Carlo 
scheme designed to model optical trap pulling experi- 
ments. Using as inputs the free energetic properties of 
DNA in relevant conformations, our model offers thermo- 
dynamic and kinetic predictions on experimentally rele- 
vant length and time scales. We operate our model in 
two 'modes', corresponding respectively to our interpre- 
tation of the force-melting and B-to-S pictures, under 
conditions designed to model experiments reported in 
Refs. [H, E3, EH ■ In 'force-melting mode' we allow each 
basepair to instantaneously adopt one of three discrete 
states: a standard helical state (B), a molten state (M) 
of two load-bearing but unhybridized strands, and an 
'unpeeled' state (U) consisting of unhybridized strands 
of which only one is load-bearing (for illustrations, see 
Figure [5]). The latter state is required to model the 
separation of one strand from another, which may oc- 
cur when breaks(' nicks') interrupt the sugar-phosphate 
backbone [1, 0, Q. Nicks are introduced in most over- 
stretching experiments as a means of relaxing torsional 
constraints. In 'B-to-S mode' we permit states B, M and 
U, and in addition assume the existence of an elongated, 
hybridized form of DNA, the S-state. We specify that 
this state is 70% longer than B-DNA, and in equilib- 
rium with the latter at 65 pN; within this framework, 
the resulting emergent dynamics of our model is largely 
insensitive to further microscopic details of this putative 
state of DNA, which we make no attempt to assert or 
infer. 

A typical experiment consists of elongating at con- 
stant rate a system consisting of A-phage DNA teth- 
ered to an optical trap. After the overstretched form 
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of DNA is observed, the elongation is reversed, allow- 
ing the molecule to recover its original contour length. 
A typical stretching-shortening cycle is reversible at low 
temperature (or high salt concentrations). As ambi- 
ent conditions are altered so as to destablise the hy- 
bridized state, force-extension data become increasingly 
hysteretic 0, [l3| . One notable characteristic of this hys- 
teresis is its asymmetry: within some temperature range 
the overstretching plateau displays little or no pulling 
rate-dependence Q [l5[ , while the behavior of the short- 
ening transition can vary strongly with pulling rate. This 
variation becomes more pronounced as temperature is in- 
creased [lH . Thus, one requirement of a viable theory of 
overstretching is the prediction of a non-hysteretic tran- 
sition at low temperature and an 'asymmetric' hysteretic 
transition at higher temperature. 

Within our model we find the following. If we as- 
sume that there exists an extended double-stranded form 
of DNA (S-DNA), and further make a minimal set of 
physically-motivated assumptions about its thermody- 
namic and mechanical properties, we find that over- 
stretching at low temperature corresponds to the in- 
terconversion of two double-stranded forms of DNA (B 
and S). Because this process involves only local free en- 
ergy barriers it is relatively rapid, and no hysteresis is 
observed. At higher temperature unpeeling can occur, 
starting from nicks in the phosphate backbone. Because 
of the interplay of the three phases B, S and U, and 
their distinct mechanical responses, we observe asymmet- 
ric hysteresis when unpeeling is considerable. We there- 
fore observe an asymmetric hysteresis that becomes more 
prominent as temperature is increased. Under conditions 
designed to model experiments reported in Refs. [l4l.ll5j|. 
we observe, following the B-to-S transition, a nonequi- 
librium high-force unpeeling. The force at which this 
unpeeling occurs is strongly rate-dependent, and agrees 
with values obtained in experiment. 

In force-melting mode our model reproduces several, 
principally thermodynamic trends observed in experi- 
ment, such as the lowering of the overstretching plateau 
with increasing temperature [3, S 0|- Detailed ki- 
netic predictions, however, do not agree with experi- 
mental findings. Unpeeling in our model within a force- 
melting framework induces a 'symmetric' hysteresis, with 
both shortening and stretching transitions falling out of 
equilibrium. The height of the overstretching plateau 
is consequently pulling rate-dependent at all tempera- 
tures. The symmetry of the hysteresis derives from the 
similarity in mechanism for strand detachment and re- 
annealing: bases open sequentially when strands detach, 
and close sequentially when strands re-anneal. These 
processes do not differ sufficiently to render the resulting 
hysteresis asymmetric. In the presence of the S-phase, by 
contrast, an asymmetry appears in force-extension data 
courtesy of the interplay of three mechanically distinct 
phases (B,U,S). However, we do not discount the possi- 
bility that a more thorough treatment of force-melting 
theory could lead, within our model, to better agree- 



ment with experiment. We summarise our findings in 
Section EH 

We draw considerable inspiration from the insightful 
theoretical treatment of overstretching of Cocco et al. [8( , 
although, as discussed later, key features of our model 
and conclusions differ from those presented in that work. 

Although most overstretching experiments are per- 
formed upon nicked DNA, and therefore permit unpeel- 
ing, our model allows us to make predictions for the case 
in which unpeeling is suppressed but the DNA remains 
torsionally unconstrained (for example, un-nicked DNA 
either in vivo or in vitro with an appropriate tethering 
arrangement). In this case, molten DNA is the only pos- 
sible unhybridized form. We find that the nonextensive 
scaling of the entropy of loops of ssDNA drives an asym- 
metric hysteresis whose qualitative character depends 
strongly on the loop exponent c of ssDNA (see Equa- 
tion (TT])). For large values of c (~ 2), the size-dependent 
free energy barrier associated with a molten loop slows re- 
annealing sufficiently to generate considerable hysteresis 
during shortening at typical experimental pulling rates. 
This hysteresis is a very sensitive function of c: for val- 
ues of c ~ 1 , we observe little hysteresis associated with 
melting at typical experimental pulling rates. The pre- 
dictions of our model for the kinetics of force-induced 
melting (when unpeeling is suppressed) therefore depend 
on how c varies with force and temperature. This de- 
pendence is currently not known. We speculate that if c 
increases with temperature (as suggested in Ref. [Hj]), we 
might again observe increasing hysteresis with increasing 
temperature. 

This paper is organized as follows. In the following 
section we discuss the possible origins of the hysteresis 
that arises during overstretching experiments. In Sec- 
tions IIIII and IIVI we introduce a simple model for over- 
stretching. We present results obtained from our model 
in force- melting mode in Sections IV Bl and IVI Al and dis- 
cuss its operation in B-to-S mode in Sections IV Al and 
IVIBI In Section IVIII we conclude by suggesting how a 
more thorough accounting of force-melting theory might 
lead within our model to better agreement with experi- 
ment, and by proposing an experiment that may permit 
the determination of the loop entropy parameter for DNA 
under tension. 



II. ENTROPY-DRIVEN HYSTERESIS 

The most striking kinetic effect observed in optical 
trap pulling experiments is hysteresis at high tempera- 
ture (e.g. above 30°C at 500 mM NaCl [l|) in force- 
extension data for the stretching and subsequent short- 
ening of a pair of initially hybridized DNA molecules. 
For a large range of temperature this hysteresis is asym- 
metric: the stretching transition appears to be similar 
at different temperatures [l3[ and pulling rates [3, EH] , 
whereas the shortening transition can exhibit pronounced 
temperature- and pulling rate-dependence. 
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Hysteresis signals physical relaxation times in excess of 
the timescale of an experiment. Common origins of such 
slow dynamics include large energy barriers due to strong 
local interactions, large free energy barriers associated 
with collective reorganization, and the intrinsically slow 
evolution of long-wavelength motion. The local interac- 
tions disrupted by DNA overstretching are noncovalent 
in character and far too weak to explain the discrepancy 
between basic rates of basepair opening (~ 10~ 8 — 10~ 4 
s [1, [13, El|) and apparent relaxation rates (seconds to 
minutes). The implicated importance of collective dy- 
namics, however, is highly unusual for systems macro- 
scopic in only one dimension, because in one dimension 
the interfaces between domains do not grow with increas- 
ing domain size. Long-ranged correlations are difficult to 
accumulate if the penalty for destroying them is micro- 
scopic, as is the case, for example, when melting a single 
turn of the double helix. Possible explanations for the 
hysteresis observed in DNA overstretching are thus few 
in number and special in character. 

In this paper we consider two such explanations. The 
dominant mechanism through which hysteresis arises in 
experiments performed to date is most likely the 'un- 
peeling' of one strand from the other, beginning from 
nicks in the phosphate backbone. This has been sug- 
gested by other authors @, 0, [1]; here, we investi- 
gate this mechanism in a framework that incorporates 
thermally-driven dynamics occurring on a broad range of 
timescales. Within our model we indeed observe hystere- 
sis when we allow unpeeling in the presence of nicks. This 
hysteresis arises from the slow, long-wavelength progres- 
sion of the unpeeled 'front' as it advances (or retreats), 
basepair by basepair, over a macroscopic distance. In the 
presence of the S-phase this hysteresis is asymmetric. For 
moderate temperatures unpeeling begins at high forces, 
at or above the force of the B-to-S transition. The slow 
conversion of S- to U-DNA at high forces is not read- 
ily seen in force-extension plots, because the mechanical 
responses of these phases in the high-force regime are 
so similar (see Section HIT]) . By contrast, the slow con- 
version of U-DNA to B-DNA upon relaxation shows a 
strong force-extension signature, because the mechani- 
cal responses of these states in the lower-force regime 
are very different. We present results obtained using the 
sequence of A-phage DNA, but sequence heterogeneity 
is not required to generate such hysteresis within our 
model. In the absence of the S-state (with or without 
internal melting), unpceling-generated hysteresis is sym- 
metric, because the mechanism of the advance and re- 
treat of the unpeeled front is the same for stretching as 
for shortening. 

Because molten DNA is thermodynamically unstable 
with respect to unpeeled DNA [8| , in the presence of nicks 
unpeeled DNA predominates over molten DNA under all 
conditions. We present simulation results corresponding 
to this case, relevant to most experiments, in Section fVl 
However, our model allows us to make predictions for 
overstretching in the absence of nicks, where the only pos- 



sible form of unhybridized DNA is the molten state. This 
scenario is of considerable interest in biological contexts. 
We find that when melting is substantial, there exists a 
possible alternative mechanism for generating hysteresis, 
originating in the conformational entropy of 'bubbles' of 
molten DNA. This conformational entropy induces long- 
ranged correlations that may persist for times compara- 
ble to the shortening protocol. For the remainder of this 
section we shall outline the mechanism by which melt- 
ing can generate hysteresis: we explore the effect of this 
mechanism upon overstretching in Section PVTl 

The entropy associated with molten regions of tension- 
free DNA does not grow linearly with the number of con- 
secutive unhybridized bases. A 'bubble' of n unpaired 
bases can be regarded as a loop of single-stranded DNA 
of length 2n, whose entropy has a contribution [2(| 
that grows sub-linearly with n: 



A5 , i 00 p(n) = -fc B cln(27 



(1) 



This nonextensive correction reflects the loss of entropy 
induced by constraining a freely-fluctuating strand to end 
at its origin (thereby forming a loop). Pinching a bubble 
in its middle incurs a statistical penalty ~ n~ c (see be- 
low). A positive loop exponent c therefore favours large 
bubbles over multiple smaller ones. It has been known 
for decades that this effect leads to a 'sharpening' of the 
melting transition [l9l |. 

The nonextensive entropy of ssDNA also has a direct 
consequence for the kinetics of melting and rehybridiza- 
tion. Assuming bulk phase coexistence between helical 
and ssDNA, the probability of nucleating a single base- 
pair of ssDNA within the helix is proportional to 
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where uq ~ 10~ 5 is the local penalty for boundaries 
between helical and molten DNA [22|. Conversely, the 
probability of nucleating a single basepair of helical DNA 
in the centre of a loop of n basepairs is proportional to 
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(3) 



From Equation we see that the free energy penalty 
for the nucleation of ssDNA is purely local, being insen- 
sitive to the size of the helix within which it attempts to 
nucleate. By contrast, from Equation ([3]) we see that the 
free energy barrier for the nucleation of double-stranded 
(dsDNA) within a region of ssDNA, which we call loop 
'pinching', grows with the size of the loop in question. 
Pinching large bubbles is very difficult. Since the clos- 
ing of a bubble from its ends is a slow process (for the 
same reason that unpeeling is slow), large molten bub- 
bles can resist complete closure for long times even when 
thermodynamically unstable. The degree of resistance to 
pinching depends strongly on the value of the loop expo- 
nent c. Estimates for c under zero tension range from 
- 1.8 120] to - 2.1 Efll. 
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We can illustrate the difference in mechanism between 
melting and reanncaling using the simple system shown 
in Figure Q] Here we show a Glauber Monte Carlo sim- 
ulation of a one-dimensional caricature of a 500 basepair 
fragment of DNA. The model consists of a blue ('hy- 
bridized') phase having energy per unit length e < 0, and 
a red ('molten') phase having entropy per unit length 
TS m > (T is temperature). In addition, the molten 
phase possesses a nonextensive entropy —cT \n(2N), N 
being the number of consecutive molten basepairs. For 
simplicity we ignore sequence heterogeneity. Starting 
from an initially hybridized state, we subject this system 
to a heating-cooling cycle. As temperature is increased, 
ssDNA becomes favoured thermodynamically, and begins 
to nucleate within the helix. At sufficiently high tem- 
perature, when most of the system is unhybridized, we 
reverse the heating protocol. As temperature falls, hy- 
bridized DNA becomes increasingly favoured thermody- 
namically, and sections of the system attempt to anneal. 
For small values of c (panel a) the free energy barrier 
to loop pinching is small; molten bubbles can close by 
pinching. In this regime melting-hybridization dynamics 
is roughly symmetric. For larger values of c (panel b) 
large loops resist pinching, and close in general slowly 
from their from their ends. Provided the heating-cooling 
cycle is sufficiently rapid, melting-hybridization dynam- 
ics is asymmetric. Dynamical symmetry is of course re- 
stored at equilibrium (panel c) , where the critical nucleus 
for melting spans the entire length of the molecule. 

This melting-hybridization asymmetry has conse- 
quences for the kinetics of conformational changes of 
DNA whenever melting is appreciable. In the context 
of DNA overstretching, we expect melting to be substan- 
tial only when unpeeling is suppressed (for example, in 
the absence of nicks). If this is so, we find that melting 
likely plays a role in both pictures we have described for 
the overstretching transition. In the force-melting picture 
the elongated form of DNA is unhybridized. Within the 
B-to-S picture, at hight temperature melting pre-empts 
the B-to-S transition. If we adopt the zero-tension form 
of the loop factor, Equation ([T]), and take c ~ 2, we 
indeed find in simulations of overstretching that when 
melting is appreciable, considerable hysteresis is gener- 
ated during the shortening transition at typical loading 
rates, regardless of whether we assume the existence of 
S-DNA. In order to fully understand the kinetic predic- 
tions of each picture of overstretching (in the absence of 
nicks), it is necessary to understand the behavior under 
tension of the the nonextensive entropy of ssDNA. 

There is currently no experimental data describing the 
force dependence of the nonextensive entropy of ssDNA. 
Calculations that neglect the excluded volume of single 
strands (see Appendix A) suggest that in the presence of 
a force / the nonextensive loop entropy takes the form 

A5ioop(n,/) = -Afe cln(2n) + k B \ng(f). (4) 

Although this entropy as a whole is tension depen- 
dent, the exponent c is not. The correction factor g(f) 




FIG. 1: Illustration of the asymmetry of melting- 
hybridization dynamics. We show microscopic configurations 
(vertical) as a function of time (horizontal) for a 500-basepair 
toy model of DNA having two phases, hybridized (blue) and 
molten (red). The sequence is homogeneous. Top panel: 
heating-cooling cycle for c = 1 (a) and c = 2 (b). We increase 
the temperature linearly with time until well above coexis- 
tence, and then reverse the protocol until the original tem- 
perature is recovered. Molten bubbles grow via nucleation, 
and close via nucleation when the parameter c is small (left 
panel). When c is large (right panel), cooling is sufficiently 
rapid that bubbles are constrained to close via domain wall 
drift. This asymmetry can give rise to melting-driven asym- 
metric hysteresis in pulling simulations, in which the vary- 
ing control parameter is tension rather than temperature (see 
Section IVlj) . Bottom panel (c): equilibrium dynamics for 
c = 2 at a fixed temperature near coexistence. With no ex- 
ternal driving, dynamical symmetry is restored. 



does not depend on n and can be absorbed into the 
local surface tension. The penalty for loop-pinching, 
Wpinch(/) = o"o 9{f) n ~ c , relative to loop nucleation, 
w nU cicatc(/) = 0o g(f) 2~ c , is essentially unchanged by 
force. Note that at forces typical of overstretching, 65 
pN, c/(65pN) ~ 10 3 , resulting in a renormalized bound- 
ary parameter of ay(65pN) = er g(65pN) ~ 10~ 2 and a 
reduced 'cooperativity' of melting. 



Calculation [23| suggests that when the excluded vol- 
ume of ssDNA is accounted for, the correction g(f) be- 
comes n-dependent. Intuition further suggests that c 
may be reduced by tension, to the extent that loops in- 
teract less strongly with a taut backbone than with a 
floppy one [24]. In simulations we vary c to identify the 
regime at which hysteresis sets in. We find that a value of 
c <~ 1 would yield a largely non-hysteretic force-induced 
melting transition at experimental loading rates. 
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III. MODEL THERMODYNAMICS 



A. Modeling force-induced melting 



To construct a model of overstretching we imagine 
coarse-graining a DNA molecule into sites spanning the 
distance along the backbone between adjacent basepairs, 
'blurring out' details on smaller scales. We illustrate this 
idea in Figure [2j Each site ('basepair') is fixed to be of 
type AT or CG (we assume fully complementary strands), 
and at any instant adopts one of the states B, M, or U. 
In the following subsection we shall introduce a fourth 
state, S. Such discrete models of biopolymers have a long 
history, beginning with the two-state model of Bragg and 
Zimm [22j, used in the 1950s to study the melting tran- 
sition of homopolymers. More elaborate discrete models 
of conformational transitions in DNA have since been 
considered, including the model of SantaLucia (25j, and 
a static three-state model of overstretching These 
models allow one to make predictions on time and length 
scales (seconds, thousands of basepairs) greatly in excess 
of those available to molecular dynamics techniques (mi- 
croseconds, tens of basepairs) [26|, albeit at a much lower 
level of detail. 

To model the fragments of A-phage DNA used in ex- 
periments we study a system of size N = 4.1 x 10 4 [l3| 
or N = 4.5 x 10 3 [H, [3 as appropriate, with sites fixed 
to be of type AT or CG according to the sequence of A- 
phage DNA [U • We denote site location by the variable 
i € [1, AT], site type by the variable <E {AT,CG}, and 
the state of a given site by the variable a <E {B,M,U,S}. 
We shall use the terms 'site' and 'basepair' interchange- 
ably. 

The equilibrium probability that a basepair adopts a 
given state a in our model is set by the free energetic 
properties of DNA in the relevant conformation. We shall 
quote energies in units of fceT at 31 OK, which we write 
as fc B T 3 i . Recall that fc B T 3 i « 0.617 kcal/mol « 4.265 
pN nm. The dynamics of our model, discussed in the 
following section, follows naturally from these statistical 
weights via a standard Glauber Monte Carlo algorithm. 

In the remainder of this subsection we shall focus on 
the force-melting theory of DNA overstretching, and al- 
low only states B, M and U. In Section IIII Bl we shall 
model the putative S-state of DNA. 

We call the logarithm of the statistical weight in our 
model, TL — — k-e,T\nW, the 'Hamiltonian': it is a free 
energy in which we imagine atomistic degrees of freedom 
have been integrated out. We separate the Hamiltonian 
into local and nonlocal pieces, H — J2^ =1 (T~t\ + ') • 
The local piece contains the extensive free energy of each 
basepair, as well as interactions between adjacent base- 



pairs. For a given site i, 

K\= - [e Tl +w B (f)]Sf 

- [TS m (T)+w M (f)} Sf 1 

- [TS m (T)+w v (f)} 6? 

+ ^BM J2 ( 5 ) 

7=M,U 

We have introduced the binary variables 8f , which take 
the value 1 when site i assumes the state a, and oth- 
erwise. The variable / is the molecular tension. 

The nonlocal piece accounts for the nonextensive en- 
tropy of molten loops, and for a loop beginning at site i 
reads 

N-i n-1 

Uf = -T £ AS loop (n, /) S^Mn II Cs- (6) 

n=l s=0 

The complementary binary variables S f = 1 — 5f take 
the value when site i assumes the state a, and 1 other- 
wise. The nonextensive loop entropy penalty can be cal- 
culated, for example, within the freely-joined chain model 
(ignoring volume exclusion), as [23j —TASi oop (n,f) = 
k-gT cln(2rt) — /cbT In <?(/). The force-dependent func- 
tion g(f) provides an approximate means of accounting 
for the change in loop entropy with force, and reads [231 ] 

^^(/,)(« ) + f |-(/,))'" m 

Here f x = fP ss /(kBT) is a reduced force, P ss « 0.7 
nm is the persistence length of ssDNA, and C(x) = 
1 - [coth(a;) - l/xf. 

In what follows we shall define the parameters of Equa- 
tion ([5]). The first three lines model the extensive free 
energies of homogeneous B, M and U phases, respec- 
tively. The fourth line describes the unfavourable free 
energy associated with interfaces between helical and ss- 
DNA . The parameters of the Hamiltonian associated 
with each phase derive from molecular free energies of two 
types: force-independent and force-dependent. We shall 
describe the force- independent parameters first. 

Force-independent free energies. B-DNA possesses 
both basepairing and stacking interactions, which we de- 
note by e n ; the variable Tj 6 {AT, CG} specifies the type 
of basepair i . We estimate these interactions from ther- 
mal melting studies. We assume that these energies de- 
pend only upon the type Ti of basepair i and not on the 
type of the neighbouring sites i ± 1. This frequently- 
made approximation can be regarded as a 'renormal- 
ization' of the neighbour-dependent stacking energies of 
more detailed models 25], by dividing into two sets in- 
teractions containing 1) any A or T bases, or 2) any C 
or G bases. This level of detail is sufficient to capture 
the essence of sequence heterogeneity, namely the extra 
stability against de-hybridization of CG over AT base- 
pairs, which in free energy terms is roughly 2 /chI^io. It 
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FIG. 2: Illustration of our model. The model is one dimensional, and consists of JV 'sites' corresponding to coarse-grained 
basepairs. Each site is of type CG or AT according to the sequence of A-DNA (sequence-dependence not shown) and at any 
instant will correspond to one of states B, M, U or S. We show an example configuration (top) of a fragment of the model, 
illustrating the conformation of DNA that we imagine corresponds to a given configuration of sites (sequence of letters). Shown 
below the configuration are schematic pictures of each state, together with their extensive and nonextensive free energies [see 
Equations 1(5} , (|6} and (|13[) ]. On the right we show free energies associated with interfaces between different states. We operate 
our model either in force-melting mode (in which we permit states B, M and U) or B-to-S mode, in which case we permit in 
addition the hypothetical S-state of DNA. 



also allows us to capture the variation with AT content 
of the melting temperature of DNA. We fix the values 
of the ener gies ecG and eat using values from the liter- 
ature 0, HHH (we assume that the entropy of melting 
liberated per basepair is S m = 12.5 /cb; see below). We 
take 

e AT = 13.17+0.19 ln(Af/0.150) 
e C G = 15.28 + 0.19 ln(M/0.150), (8) 

in units of /cbT3io- Here M is the molar NaCl concentra- 
tion (we quote the salt-dependent correction from 0). 
We choose as a salt concentration either 150 mM NaCl 
(so as to model the experiments of Refs. [H], HH) or 500 
mM NaCl (in order to model experiments reported in 
Ref. 0). 

In addition, in B-DNA there exists an interaction be- 
tween each nucleotide and a nucleotide on the oppos- 
ing strand four basepairs distant. This interaction is 
lost upon disruption of the helix [22|, H3, We ac- 

count for this effect by assigning an interaction parame- 
ter £bm ~ 11-51 k^T [1, H3] to neighbouring B- and non- 
B sites (fourth line of Equation ([5])). This interaction 
corresponds to a junction penalty, familiar from Bragg- 
Zimm 22] models, of er = cxp ( — £bm/ ikuT)) — 10~ 5 . 

The statistical weight of the M-phase receives a contri- 
bution from the extensive entropy per basepair (bp) lib- 
erated upon melting (under zero tension), 5 m (T), which 



we take to be 0, |1] 

S m (T) = S m (T m ) + i^Pj < W 

where S m (T m ) = 12.5 /cb- The second term (in which 
we have absorbed both entropic and enthalpic contribu- 
tions Q) derives from the change in heat capacity upon 
melting, quoted as AC P « 60 ± 10 cal/(mol-K-bp) 
We set AC P = 60 cal/(mol-K-bp). The specific heat 
correction to the free energy increases from zero at 
the melting point to about lfcsT3io at 0°C. The ex- 
perimental melting temperature T m for A-phage DNA 
(which has ^50% CG content) has been fit by 0, Q 
T m (°C) w 104.4 + 16.89 \og 10 M. This formula gives a 
melting temperature of about 90°C at 150 mM, in rea- 
sonable agreement with the zero-force melting point of 
our model (about 87°C for c = 2; data not shown). 

We assume the unpeeled state (U) to be unhybrizided 
DNA consisting of one load-bearing single strand; the 
other strand is free to fluctuate [|[. We take its exten- 
sive entropy to be that of M-DNA, and we assume that 
U-DNA has no nonextensive entropic correction. In our 
model, any unhybridized site connected to at least one 
nicked site by an unbroken sequence of unhybridized sites 
is defined to be unpeeled; if not so connected, it is defined 
to be molten. We assume that there is no energetic bar- 
rier between the M and U states Q. We set the energy 
of U-B junctions equal to those of M-B junctions. 

Force- dependent free energies. The terms w a (f) in the 
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first three lines of Equation ([5]) describe the elastic free 
energies of the B, M and U states of DNA. One can de- 
termine these energies by integrating, as a function of 
pulling force, the experimentally-measured extension per 
basepair x a (f) of each state i.e. 



df'x a (f) 



(10) 



For forces of interest, the extension per basepair of B- 
DNA as a function of force / is give n ap proximately by 
an extensible worm-like chain model |2£ 



xn(f) = x B 1 



1 ( kBT 

2 \fPs 




(11) 



with contour length xb — 0.34 nm/bp, persistence length 
Pb = 50 nm, and stiffness Sb = 1300 pN. We bound the 
above expression (which breaks down at very small forces 
/ < 0.1 pN) from below by zero. 

We take the mechanics of M-DNA to be that of two 
pieces of load-bearing ssDNA, for which we adopt the 
parameterization of Ref. Q , 



,(/) 



ailog(///i 



1 + a3e 



-f/h 



a 2 ~ — 



(12) 



Here x ss = 0.34 nm, a\ — 0.21, a 2 — 0.34, 03 = 
2.973 + 0.5129 InM, h = 3.7 x 10~ 3 pN, / 2 = 2.9 pN, 
and fs — 8000 pN. This formula accounts for the increase 
in persistence length of ssDNA under low salt conditions. 
We bound this expression from below by zero. We find 
that the extensible freely-joined chain fit to ssDNA [30| . 
which does not account for the increase in its persistence 
length at low salt concentrations, leads to numerically 
different melting forces at nonzero tension, but does not 
affect the qualitative predictions of our model. We ne- 
glect the effects of secondary structure, which for load- 
bearing strands we expect to be small. 

We assume that M-DNA consists of two noninter- 
acting, load-bearing strands, and so we take im(/) = 
x ss {f/2) and u>m(/) = 2w ss (//2). We asssume that U- 
DNA consists of a single load-bearing strand [8J , and so 
take iu(/) = x ss (f) and w v (f) = w ss (f). 

In principle there is a force-dependent contribution to 
the entropy of ssDNA arising from the 'pulling out' at 
high forces of its backbone degrees of freedom. This effect 
is negligibly small, however, as may be shown within the 
freely-jointed chain model of ssDNA. 

Finally, we approximate the nonextensive loop entropy 
of ssDNA via Equations © and Q [23]. This approx- 
imation neglects the excluded volume of single strands. 
In view of this, we regard the loop exponent c as an un- 
known and vary it between and 2.1. 

We shall examine the predictions of our model both in 
the presence and the absence of nicks. In the presence of 
nicks (the usual experimental situation) we permit un- 
peeling. Within our model we find that molten DNA is 
always thermodynamically unstable to unpeeling. This is 
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FIG. 3: Force (vertical axis) at which M- (red) and U-DNA 
(blue) are in equilibrium with B-DNA, as a function of tem- 
perature (horizontal axis). This calculation is based on the 
extensive free energies of the Hamiltonian, given in Equa- 
tion The lower pair of bold curves refers to AT sites; the 
upper pair refers to CG sites. The thin line was computed 
in the absence of the specific heat correction shown in Equa- 
tion This correction markedly lowers the overstretching 
force at low temperature, giving much better agreement with 
experiment. With no free energetic barrier between M and 
U, molten DNA is always thermodynamically unstable to un- 
peeling [8j. In the presence of nicks (the usual experimental 
situation) we therefore find unpeeled (U) DNA to be the pre- 
dominant form of unhybridized DNA. 



so by virtue of the unpeeled state being longer at a given 
force than the molten state (a single load-bearing strand 
extends more than do two strands each bearing half the 
load) . A similar conclusion was reached in Ref. [8j . The 
free energy of U-DNA can fall below that of the B-state 
at a lower force than can the free energy of M-DNA. 
There is therefore a range of T and / for which unpeel- 
ing can cross CG regions (and therefore proliferate over 
long distances), while melting cannot. 

We demonstrate the instability of melting to unpeeling 
in Figure [3j in which we show the force (as a function of 
temperature) at which B-DNA is unstable to unpeeling 
and to melting. For nicked molecules the predominant 
form of unhybridized DNA is therefore unpeeled. We 
shall examine the predictions of our model in this case 
in Sections IV Al (B-to-S mode) and IV Bl (force-melting 
mode). 



B. Modeling S-DNA 

To model the picture of overstretching as a transition 
to an elongated form of dsDNA, we extend our model to 
include a fourth state, S-DNA @, 0, ©, M, HI G3- S- 
DNA, whose existence is debated, is naturally less well- 
characterized than either helical or molten DNA. We 
must therefore make several assumptions about its ther- 
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modynamic and mechanical properties. However, within 
our model a minimal set of physically- motivated assump- 
tions results in a kinetics whose qualitative features are 
robust to changes in model details. Our key assumptions 
in this section and the section following arc 1) that S- 
DNA exists; 2) that it is double-stranded (and so its free 
energy possesses no long-ranged entropic component); 3) 
that it is 70% longer than B-DNA, and 4) that it is in 
thermal equilibrium with the latter at 65 pN. Within this 
framework the qualitative kinetic details that emerge de- 
pend very weakly on the precise numerical values of our 
S-DNA parameters. 

The B-to-S picture must be considered tentative until 
there exists a direct observation of this state of DNA. The 
force-melting picture, by contrast, is a first-principles ap- 
proach with considerable predictive power, building on 
well-known properties of helical and molten DNA. In 
this study we do not propose or defend a likely atom- 
istic model for S-DNA. More detailed computer simula- 
tions [26| are required to investigate this issue. Instead 
we ask what are the likely kinetic consequences of assum- 
ing that overstretching can involve two double-stranded 
forms of DNA. This assumption has been made by sev- 
eral authors 0, B S [M El El • We extend our model 
accordingly. 

In 'B-to-S' mode we permit basepairs to adopt the 
forms B, M and U as previously, and in addition permit 
the adoption of a fourth state, S. We define the ther- 
modynamic properties of this state via the Hamiltonian 
ffs=EL«i, where 

nf = -[X n e n +w s (f)]Sf 

7=M,U 

+ ^bs (SfSf +1 + SfSf +1 ) . (13) 

The first line of Equation (TT3"|) describes the extensive 
free energy of S-DNA. The second term in the first line 
of Equation (fT3|) . ws(f), is the elastic free energy of S- 
DNA. Based on the idea that S-DNA is double stranded 
and stiff [8j, we assume that it can be parameterized by 
the wormlike chain in a similar fashion to B-DNA. We as- 
sume that a?s(/) has a form similar to Equation (fTT|) . with 
subscripts 'B' replaced by subscripts 'S'. We assume a 
contour length xa = 0.58 nm/bp, and a stiffness constant 
S§ = 2730 pN [8|. Lacking experimental measurements 
of S-DNA flexibility, we assume its persistence length to 
be equal to the persistence length of B-DNA. The emer- 
gent dynamics of our model is qualitatively insensitive to 
variation of the latter two parameters, which determine 
the shape (but not the overall scale) of the equilibrium 
extension curve for S-DNA. 

The factor X Ti , where Ti £ {AT, CG}, quantifies the di- 
minished base-pairing-stacking energy of S-DNA relative 
to B-DNA. We determine this loss by requiring phase co- 
existence between the B and S phases at a pulling force 



of f = 65 pN, i.e. we fix X Ti by setting 

A Ti e Ti + w s (f ) = £ Ti + w B {fo)- (14) 

We obtain Aat ~ 0.74 and A C g ~ 0.77. We thereby 
assume that the B-to-S coexistence forces for AT and 
CG sections are equal. It is possible that this equality 
does not hold for real DNA 0, 0| , but conclusive data do 
not yet exist. Indeed, there is evidence suggesting that A- 
DNA (50:50 AT:CG) and CG-rich DNA both overstretch 
at 65 pN [l5|. The qualitative kinetic features we shall 
describe are robust to variations in these forces of tens of 
piconewtons. 

The second line of Equation (TIB"]) describes the ener- 
getic penalty for interfaces between S-DNA and molten 
DNA. We assume this to be identical to the B-M junction 
penalty (and to the B-U and S-U penalties). 

The third line of Equation (|13[) describes the ener- 
getic penalty for interfaces between S-DNA and B-DNA. 
The parameter ess controls the slope ('cooperativity') 
of the overstretching (B-to-S) plateau at low tempera- 
ture. In Ref. [|| this was set to 3/cbT 1 - We found that 
this value, in concert with our wormlike chain parame- 
terization for S-DNA, gave an overstretching transition 
noticably less cooperative than that observed in experi- 
ment (in Ref. [|[ a different elastic behavior was assumed 
for S-DNA). We choose instead to set the junction energy 
to ees — 5fcBT3io, in order to give an overstretching co- 
operativity similar to that observed in experiment. This 
choice is without theoretical foundation, and should be 
regarded as a fit to experimental data. However, the ki- 
netic features we observe are insensitive to variation of 
this parameter over a range of tens of U-qT] the role of 
ees is to set the slope of the B-to-S plateau. 

IV. MODEL DYNAMICS 

With the thermodynamics of our model determined, 
we investigate its predictions for the dynamics of over- 
stretching. We use a Monte Carlo (MC) dynamics de- 
signed to mimic pulling at a constant speed, appropriate 
for experiments carried out using optical tweezers. 

We consider an N basepair stretch of model A-phage 
DNA. We choose N to be either 4.1 x 10 4 (when 
modeling the experiments of Ref. [l3j|) or 4.5 x 10 3 
(when modeling experiments presented in Refs. [H], EH). 
We describe its state at a given time t by the vec- 
tor {ai(t), . . . , cti(t), . . . , aAr(i)}, where ai{t) denotes the 
state a £ {B,M,U,S} of basepair i. To this we attach a 
Hookean spring designed to model the optical trap. We 
require that the system (molecule plus trap) be in me- 
chanical equilibrium at all times, i.e. that the tension of 
the trap and at all points on the molecule have a uniform 
value, /. We thereby assume that mechanical transduc- 
tion of force (including propagation of torque) is much 
more rapid than molecular changes of state. We write the 
molecular extension as £(f,t) = J2 a n a(t) x a (f), where 
n a (t) denotes the number of basepairs in state a at time 
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t, and x a (f) is the extension per basepair of state a as 
a function of force. The optical trap bead is moved / /k t 
from the trap centre, where fc t =0.1 pN/nm is the trap 
force constant. 

In optical trap experiments the total extension of the 
molecule and trap, L(t), is increased at constant speed, 
L(t) = v t. In our simulations we impose a total length 
of this nature, and calculate the resulting tension / and 
molecular extension £(f, t) by requiring that the imposed 
extension equal the extensions of the molecule plus trap, 
L(t) = £(f,t) + f /kt. The relative extensions of the 
molecular states are fixed by the constraint of mechanical 
stability, / = /bOb) = fv(xxj) = 2/m(^m) = fs{xs), 
where the f a (x) are found by inverting the functions 
x a (f). 

In pulling simulations we impose fixed boundary con- 
ditions, with sites 1 and N constrained to be B-form 
(modeling 'clamped' ends). We start each run from a 
fully helical configuration at imposed length Lq = 0.317V 
nm, which implies Nxsit = 0) + /bO^b)/^ = Lq. This 
relation determines the initial molecular tension f(t = 0) 
and the initial extensions x a (t = 0). We equilibrate the 
system at this fixed, small extension before beginning the 
pulling protocol. To perform a pulling simulation we in- 
crement the global length according to L(t) — vot, with 
t measured in units of Monte Carlo sweeps (on average 
one update per site). When the force exceeds 120pN (600 
pN in Figure[S]), we reverse the loading rate (vq — » — Vq). 
A Monte Carlo step consists of a proposed change in the 
state on — > a[ of a randomly chosen site i, accepted with 
the Glauber probability 

Pacc(a* eO = [1 + exp^AT^f 1 . (15) 

Here /3 = (fc B 7T\ and 

AHi = H [{an,..., a' i ,...,a N }] 

- H [{ax, . . . ,ati, . . . ,a N }} (16) 

is the change in Hamiltonian Ti upon the change in state 
of basepair i from a,i to a[. The thermodynamic param- 
eters of our model enter the rates for changes in basepair 
state via the Hamiltonian H, Equation ([5]). Recall that 
the Hamiltonian is nonlocal by virtue of the loop entropy 
term for ssDNA. Because Glauber transition rates satisfy 
detailed balance [36|, our simulated DNA molecule will 
evolve toward its equilibrium state at a given force /. 
Whether it reaches this equilibrium state depends on the 
ratio of rates for changes in molecular state and pulling: 
very rapid pulling gives the molecule little chance to equi- 
librate at a given tension. We assume that the fundamen- 
tal timescale on which individual bases attempt to change 
state ('basepair closing time') is Tf = 28/xs [l8j . chosen 
to agree with basepair fluctuation timescales measured 
using fluorescence correlation spectroscopy (T3, [HI ■ This 
is then the timescale for a single Monte Carlo sweep. We 
vary the pulling rate over the range 150 to 4000 nm/s, 
appropriate for experiments reported in Refs. [l3L [Til flql . 



The kinetic data we present are qualitatively robust to 
a change in fundamental timescale of an order of magni- 
tude. Timescales derived from FCS data associated with 
the closing of individual basepairs range from Tf ~ 10~ 4 
to Tf ~ 10~ 5 s [l7|; our results are qualitiatively un- 
changed if we adopt as our timescale any value within 
this range. We find that within our model the fundamen- 
tal timescale of Ref [18j (Tf = 28/is) gives good agreement 
with an experimentally observed rate-dependent unpeel- 
ing transition [3, [HI (see Figure [5]) . 

The question as to what constitutes the 'best' timescale 
for coarse-grained models such as the one we study here 
is an open question. Nuclear magnetic resonance (NMR) 
data imply that molecular fluctuations occur on a much 
shorter timescale, 10~ 7 to 10~ 8 s @,[l7[, than those asso- 
ciated with FCS data. However, the authors of Ref. [ItJ 
suggested that NMR fluctuation timescales might be sen- 
sitive to fluctuations smaller than those necessary to ef- 
fect the opening or closing of a single basepair. Further- 
more, a simple kinetic model employing a basepair clos- 
ing time of 28 /its was shown in Ref. [l8[ to give reasonable 
agreement with experimental FCS data. 

In Ref. [8J] a base-closing timescale of 10~ 8 s was 
adopted and used to derive, for a simple kinetic model, 
predictions for the nature of the hysteresis seen dur- 
ing DNA overstretching. The model concerned the ad- 
vance and retreat during overstretching of the unpeeled 
'front' (the junction between unpeeled and hybridized 
DNA) over an inhomogencous sequence. The authors 
showed that considerable 'roughness' in force-extension 
data arose during shortening as a consequence of the 
'stick-slip' motion of the front as it traversed the inho- 
mogeneous sequence. Although the timescale chosen in 
that work is seemingly at odds with our choice of a much 
slower FCS-derived timescale, we believe the apparent in- 
consistency can be explained by the difference in dynam- 
ical procedures used in our paper and in Ref. § . In this 
paper we use a heat bath (Glauber) dynamics in which 
thermal fluctuations permit the system to explore its free 
energy landscape. As a consequence, the emergent dy- 
namics of our model incorporates processes that occur on 
a broad range of timescales, including those that require 
the crossing of large free energy barriers (called 'acti- 
vated processes'). By contrast, the dissipative dynam- 
ics employed in Ref. Q permits no processes involving 
thermal fluctuations, instead corresponding to a mono- 
tonic descent to the local free energy minimum. Ther- 
mal processes are considered in an approximate fashion 
by smoothing (averaging) the local sequence. In the con- 
text of the model discussed in [8j , this dynamics amounts 
(for a given global extension) to relaxation of the un- 
peeled front toward the local free energy minimum of the 
sequence. As total extension is incremented, this local 
minimum is changed, and the unpeeled front may re- 
lax again to the new local minimum. We demonstrate 
in Appendix B that this dissipative dynamics is largely 
insensitive to the fundamental timescale employed, pro- 
vided that the unpeeled front has time to reach the lo- 
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cal free energy minimum. We show that this is so for 
a given set of conditions for timescales ranging from 
T f = 2.8 x 10~ 5 s to r f = 10~ 8 s. We therefore believe 
there to be no direct contradiction between our chosen 
Monte Carlo timescale of 28 /is and the 10 -8 s adopted 
in Ref. [|[ as the timescale for a dissipative dynamics. 
The difference in fundamental timescales is not evident 
in force-extension data because of the different dynami- 
cal protocols used. In the context of our model, we find 
the FCS timescales to give much better agreement with 
experiment when considering the slow, high-force unpeel- 
ing transition following the initial overstretching plateau 
(Refs. HUll and Figure EJ). 

In parallel with proposed changes of molecular state we 
propose N times per MC sweep a change 5f in the global 
tension, accepted if the quantity \ ^2 a n a x a + f /k t — L(t)\ 
does not increase. This frequent adjustment of force en- 
forces the condition of mechanical equilibrium. 

The behavior of U-DNA is governed by the following 
rules. We allow unpeeling to start only from nicks, which 
we place randomly in every simulation (in some simula- 
tions we consider un- nicked DNA). We consider all nicks 
to be on the same strand, ensuring that one strand is 
always load-bearing. A nick is made between randomly- 
chosen neighbouring sites, i n and i n + 1. Unpeeling can 
proceed left from i n and right from i n + 1. Any unhy- 
bridized site connected to at least one nicked site by an 
unbroken chain of unhybridized sites is defined to be un- 
pcclcd; if not so connected, it is defined to be molten. 
The effect of local changes of state associated with a 
change in state of a single site may be non-local. For 
example, annealing an unpeeled strand of length l\ at a 
position £2 < £1 from its end changes a region of length 
£1 — £ 2 from unpeeled to molten. Such a process is gen- 
erally very unfavourable, however, being suppressed by 
the thermodynamic predominance of U- over M-DNA. 
If unpeeled and molten regions merge, the latter is sub- 
sumed by the former (as is appropriate given our defini- 
tions of 'unpeeled' and 'molten'). If an unpeeled front 
advances from its starting nick to an adjacent nick, the 
corresponding unpeeled region is considered to be com- 
pletely detached, and is not permitted to change state 
for the remainder of the simulation. 

Kinetic data shown in the following sections are repre- 
sentative individual trajectories, selected at random from 
a large number of such trajectories. On some plots we 
show symbols denoting equilibrium force-extension val- 
ues. Each equilibrium value was calculated by averaging 
four independent simulations run at zero loading rate; 
error bars are smaller than symbols. In Figures [6] and [8] 
we show in addition equilibrium force-extension data for 
nicked systems in force-melting mode. Because in this 
case complete enumeration of microstates is straightfor- 
ward, these data were calculated by evaluating numeri- 
cally the partition function using the Hamiltonian of our 
model and the sequence of A-DNA, assuming a nick be- 
tween basepairs 2 and 3. 



V. SIMULATION RESULTS IN THE PRESENCE 
OF NICKS 

A. B-to-S mode 

We observe a rich dynamics in B-to-S mode, due to the 
interplay of two hybridized states (B and S) and, chiefly, 
one unhybridized state (U) . In Figure [4] (left panel) we 
show the temperature dependence of the overstretching 
behavior of once-nicked A-DNA at 500 mM NaCl and 
pulling rate 1000 nm/s, chosen to model the experiments 
presented in Ref [1J|. We set the loop exponent c to 
2.1 [33]. At low temperature, well below the zero-tension 
melting point of our model at this salt concentration 
(T m « 87°C), the hybridized state is stable with respect 
to the unhybridized state over the full range of force. As 
determined by the thermodynamics we have imposed, we 
observe a transition from a B-rich to an S-rich phase at 
the B-to-S overstretching force of 65 pN. Because the in- 
terconversion of two hybridized states involves only local, 
and therefore modest, free energy barriers, we observe a 
reversible dynamics associated with the B to S transition. 

As temperature increases, partial unpeeling begins to 
occur. U-DNA is longer than B-DNA at tensions > 10 
pN, and so becomes increasingly favoured thermodynam- 
ically as tension increases. In this regime we observe 
the progressive increase with temperature of the degree 
of an hysteresis that is asymmetric in character: the 
overstretching force remains pulling rate-independent for 
moderate temperatures, whereas the shortening behav- 
ior is strongly rate-dependent. The degree to which we 
observe hysteresis at a given temperature is similar to 
that seen experimentally in Ref [13]. In our model, this 
hysteresis arises from the slow dynamics induced by the 
macroscopic displacement of the nicked front as it moves 
basepair-by-basepair. The origin of the hysteresis asym- 
metry is the interplay between the phases B, S and U, 
combined with the mechanical behavior of each. We clar- 
ify this point in Figure |4] (right panel) by showing for 
30° C the percentage of the molecule in unpeeled form. At 
this temperature the initial transition is B-to-S, occurring 
in equilibrium at 65 pN. We also observe a slow unpeel- 
ing (conversion of B to U and later S to U) slightly in 
arrears of the B-to-S transition. However, this unpeeling 
is largely masked in force-extension data by the similarity 
of the S and U force-extension relations in the relevant 
regime. By contrast, the slow annealing (conversion of U 
to B) observed during shortening is made evident by the 
distinct nature of the U and B force-extension relations 
at low forces. 

As temperature increases further, unpeeling causes the 
overstretching force to drop (above 30° C in Figured]). 
As hybridized DNA becomes less stable with respect to 
unhybridized DNA, force-induced unpeeling (and hence 
molecular elongation) can occur at tensions at or below 
65 pN, concurrent with or pre-empting the B-to-S tran- 
sition. In this event the plateau height falls, and be- 
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FIG. 4: B-to-S mode: Effect of temperature. Simulated force-extension data for once-nicked A-phage DNA at 500 mM NaCl, 
chosen to model experiments presented in Ref. Left panel. We show data for pulling rate 1000 nm/s at T — 10, 20, 30, 
40 and 60° C (offset left to right), together with equilibrium data for 10 and 30° C (circles). At low temperature overstretching 
is a transition from B- to S-DNA. Increasing temperature causes unpeeling, resulting in the onset of 'asymmetric' hysteresis 
and a drop in plateau height. The progressive increase with temperature of the degree of hysteresis results from a competition 
between the entropy of U-DNA and the basepairing energy of S-DNA. The overstretched state at high temperature is largely 
U-DNA. The grey lines here (and in subsequent figures) are equilibrium force-extension data for pure states. Right panel. 
For 30°C we show data obtained at pulling rates of 1000, 2000 and 4000 nm/s, together with the percentage of unpeeled DNA 
during stretching and shortening at the slowest and fastest rates (top panel; arrows and line colors correspond to those on 
force-extension curves). The slow initial unpeeling is largely disguised by the similarity of the S and U force-extension relations 
in this regime (see, however, Figure [5}. By contrast, the slow annealing is made evident by the distinct nature of the U and 
B force-extension relations at low forces. In this example increasing the pulling rate reduces the degree of hysteresis observed, 
because one strand has less time to detach from the other. 



comes rate-dependent. When considerable unpeeling is 
observed the plateau also 'roughens', fragmenting into 
successive mini-plateaux (left panel of Figure [5J data for 
60 °C). This effect reflects the slow, punctuated dynam- 
ics of the unpeeled front as it crosses an heterogeneous 
sequence. 

The degree of unpeeling-driven hysteresis in B-to-S 
mode depends not only on temperature, but also on the 
position of the nick and on the pulling rate. If the nick, 
which we place randomly in every simulation, occurs in 
a CG-rich region (the first 23 kbp of A-DNA are roughly 
15% richer in CG basepairs than in AT basepairs), then 
unpeeling faces a greater barrier to proliferation than had 
the nick been placed in an AT-rich region. Thus in Fig- 
ure 2] while hysteresis in general increases with temper- 
ature, the effect of random nick placement can perturb 
this trend (see results for 10°C and 20°C). Hysteresis is 
in addition a non-monotonic function of the pulling rate. 
At very low pulling rates the molecule can remain in equi- 
librium. At a pulling rate of 1000 nm/s, hysteresis can 
be appreciable (see result for 30° C, Figure [J). At higher 
rates (2000 and 4000 nm/s; see Figure |H right panel), 
the degree of hysteresis decreases, because unpeeling has 
less time to proliferate. Again, the effect of random nick 
placement can disrupt this trend. 



More than one transition may be observed in B-to- 
S mode. Motivated by experiments reported in Refs [Til . 
[lH , we perform pulling simulations at salt concentrations 
of 150 mM up to a maximum force of 600 pN. Because 
at high forces ssDNA is more easily extended than our 
model of S-DNA, we observe for a range of temperature 
a nonequilibrium unpeeling transition at very high forces 
(upwards of about 150 pN), following the 65 pN equi- 
librium B-to-S transition. A similar interpretation was 
given to the two-stage overstretching in Refs. 
We show representative data in Figure[5](left panel). The 
second transition (S-to-U) occurs out of equilibrium, as 
evidenced by the pulling rate-dependence of the transi- 
tion force. If the terminal pulling force is large enough 
that this transition is observed, the stretching-shortening 
cycle becomes hysteretic. This is true even at tempera- 
tures sufficiently low that the B-to-S transition shows no 
hysteresis when the pulling is reversed before the force 
exceeds ^150 pN. In Figure^] (right panel) we show that 
for fixed pulling rate the unpeeling transition can at high 
temperature pre-empt the B-to-S transition, leading to 
a drop in plateau height. An apparently similar sce- 
nario was observed experimentally in Refs. [3, EH , and 
was discussed theoretically in Ref. [8j. Note that while 
the high-force unpeeling in our model is strongly rate- 
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FIG. 5: B-to-S mode: high-force transition. Simulated force-extension data for a once-nicked 4.5 kbp fragment of A-phage DNA 
at 150 mM NaCl, chosen to model experiments reported in Refs. [l4l.[l5|. Left panel. Effect of pulling rate at 20° C. We employ 
pulling rates of 3000, 1500, 700 and 150 nm/s. Following the B-to-S transition at 65 p N, we observe a strongly rate-dependent 
unpeeling transition at forces in good agreement with data presented in Refs. [l4l . HEf . We infer from this comparison that the 
correct fundamental timescale for our coarse-grained model lies in the 10-100 /xs range characteristic of fluorescence correlation 
spectroscopy data [l7l Il8|]. We show also the fraction of unpeeled DNA for pulling rates 150 and 1500 nm/s (top panel; arrows 
and colors correspond to those of force-extension data). Right panel. Effect of temperature at pulling rate 150 nm/s. We 
show data for temperatures 0, 30 and 50° C. At sufficiently high temperature the unpeeling transition pre-empts the B-to-S 
conversion, lowering the overstretching plateau. 



dependent, the initial B-to-S transition is not, occurring 
at 65 pN for all pulling rates from 150 to 3000 nm/s. 
Similar behavior was seen in experiment [l4l [FoT ] . 

The force at which we observe the S-to-U transition 
varies strongly with pulling rate, but the values we ob- 
tain are similar to those found in experiment at the same 
pulling rates (e.g. Figure 1 of Reference [3). Assuming 
that our proposed mechanism for this transition is cor- 
rect, we infer from this comparison that the basic fluctu- 
ation timescale we have adopted (28 /zs [18|) is physically 
appropriate for our coarse-grained model. 

In these high-force simulations we assume a naive ex- 
trapolation of all elasticity data to large forces. We 
also dispense with g(f), the force-mediated correction 
to the loop entropy of ssDNA, because retaining this fac- 
tor leads to a change in sign of the surface tension of 
dsDNA-ssDNA junctions at very large forces. 



B. Force-melting mode 

In force-melting mode our model can access only states 
B, M and U, and corresponds to our interpretation of the 
force- melting theory of DNA overstretching [1, H, 0, U H] • 
In Figure[6]we show data obtained under conditions iden- 
tical to those used to generate Figure 01 chosen to model 
the experiments reported in Ref . [131 ] . As the force on the 
DNA increases, B-DNA becomes unstable to the unhy- 



bridized forms M and U. Stretching results in the com- 
plete unpeeling of the DNA. This is so because of the 
thermodynamic predominance of U- over M-DNA, and 
also partially because wherever these two phases meet the 
former subsumes the latter (as is appropriate given our 
definitions of 'molten' and 'unpeeled'). Because the dy- 
namics of the unpeeled front is slow, accompanying this 
unpeeling is considerable hysteresis. However, unlike in 
the presence of the S-state, unpeeling-driven hysteresis in 
force-melting mode is symmetric: the stretching transi- 
tion lies roughly as far from equilibrium as does the short- 
ening transition. This may be seen clearly by comparing 
kinetic and equilibrium data in Figure [51 This symme- 
try derives from the fact that strands detach or re-anneal 
in a similar fashion, namely, via the sequential advance 
or retreat of the unpeeled 'front' [33|: the kinetics of 
strand detachment is no less sluggish than that of strand 
re-annealing. The advance and retreat of the unpeeled 
front occurs by the same mechanism in the presence of 
S-DNA, but there an asymmetry in force-extension data 
is induced by the interplay of three phases, B, S and U. 

The sluggishness of the unpeeling transition confers 
upon overstretching data an hysteresis at all temper- 
atures. The symmetry of the hysteresis means that 
both shortening and stretching occur out of equilibrium, 
and therefore that the plateau height exhibits a rate- 
dependence at all temperatures. By contrast, as dis- 
cussed in the previous subsection, there is in experiment 
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FIG. 6: Force-melting mode: Effect of temperature. Simulated force-extension data for once-nicked A-phage DNA at 500 
mM NaCl, chosen to model experiments presented in Ref. [l3| ]. Left panel. We show data for pulling rate 1000 nm/s at 
T = 10, 20, 30, 40 and 50°C (offset left to right), together with equilibrium data for 10 and 30°C (circles). Overstretching 
corresponds to the conversion of B- to U-DNA, with accompanying hysteresis. This hysteresis is symmetric, in the sense that 
the stretching transition lies as far from equilibrium as does the shortening transition. There is no appreciable change with 
temperature in the degree of hysteresis, because no competition between U- and S-DNA can arise (we do not permit the latter 
state in force-melting mode). Right panel. For 30°C we show data obtained at pulling rates of 1000, 2000 and 4000 nm/s, 
together with the percentage of unpeeled DNA during stretching and shortening at 1000 nm/s (top panel; arrows and line color 
correspond to those on force-extension curve). The slowness of unpeeling confers a rate-dependence upon the overstretching 
plateau. 



a range of temperature over which the transition is re- 
versible and the overstretching plateau appears to be 
pulling rate-independent [HI, O, EBl ■ In addition, there 
exists in force-melting mode no analogue of the high- 
force transition [3, EH observed in B-to-S mode (Fig- 
ure , because the unpeeled state is the overstretched 
state: there is no additional phase to which a transition 
may be made. 

In order to prevent complete unpeeling of the DNA 
in the presence of nicks, we would require within our 
model a highly collective frcc-cncrgetic or kinetic 'barrier' 
between molten and unpeeled DNA, or a different notion 
of what are the 'molten' and 'unpeeled' states of DNA. 
We discuss this possibility in Section IVIII 



VI. SIMULATION RESULTS IN THE ABSENCE 
OF NICKS 

A. Force-melting mode 

In this section we examine the predictions of our model 
in the absence of nicks. Because we forbid fraying from 
the ends of our simulated DNA molecule (the first and 
final basepairs are constrained to remain hybridized for 
the duration of the simulation), unpeeling cannot occur. 
Nicks allow one strand to rotate freely about the other 



during stretching. In the absence of nicks, we nonetheless 
assume that the DNA remains torsionally unconstrained; 
if torsionally constrained, the overstretching transition 
occurs in excess of 100 pN [lj|. Torsional freedom in the 
absence of nicks might be engineered by using a careful 
tethering arrangement. 

When unpeeling is suppressed, the only form of un- 
hybridized DNA accessible is the molten state. Con- 
sequently, as shown in Figure [7J we observe in force- 
melting mode a transition from B-DNA to M-DNA upon 
stretching. This transition occurs at roughly 70 pN at 
room temperature, as per the thermodynamic analysis of 
Refs. @, y, Note that omitting the specific heat cor- 
rection in Equation ([9]) leads to much larger overstretch- 
ing forces at low temperature (see Figure [3]) [23] . The 
plateau would be lowered still further by adopting larger 
values of AC P still within experimental uncertainty (e.g. 
w 70 cal/(mol-K-bp)). 

When melting is appreciable we observe hysteresis 
whose extent is strongly dependent upon the value of 
the loop exponent c. In Figure [7J we show data for un- 
nicked A-DNA at fixed temperature (25° C) and salt con- 
centration (150 mM NaCl), for varying loop exponent c. 
While the apparent sharpness of the equilibrium transi- 
tion (circles) is modified only weakly by variation in the 
loop exponent (this sharpness is partially counteracted 
by sequence heterogeneity), the degree of hysteresis is 
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FIG. 7: Force-melting mode in the absence of nicks: Effect 
of varying loop exponent. Simulated force-extension curves 
for un-nicked A-phage DNA at 150 mM NaCl, T = 25° C, 
and pulling rate 1.8 /im/s. The curves correspond to loop ex- 
ponents c = 0, 1.25, 1.5 and 1.8 (offset horizontally), and we 
show equilibrium data in three cases (symbols) . Hysteresis in- 
creases as c increases, even though the equilibrium fraction of 
M-DNA decreases with increasing c. Note that the 'coopera- 
tivity' (or sharpness) of the force-melting transition decreases 
as c decreases. For the largest value of c we show the effect 
of pulling rate upon the degree of hysteresis. 



strongly dependent upon this parameter. For values of c 
appropriate at zero tension (c ~ 2) we see considerable 
hysteresis at typical experimental rates of pulling. 

In our model, this hysteresis is driven by the devel- 
opment of large molten bubbles (see Section HI)) . At 
sufficiently high tension, molten DNA nucleates and 
spreads. The nonextensive loop entropy rewards bub- 
ble coalescence. During shortening, as tension decreases, 
B-DNA becomes thermodynamically favoured over M- 
DNA. However, for large values of c bubbles resist pinch- 
ing, and the nucleation of B-DNA within M-DNA is sup- 
pressed. Large molten bubbles thus close in general from 
their ends. Because this is a gradual process, requiring 
a collective relaxation over often large distances, such 
bubbles can persist far into the shortening transition, 
preventing the molecule from reaching equilibrium. For 
smaller values of c bubble pinching is much less costly, 
and equilibrium is attained more readily [3c| . 

We find that melting-driven hysteresis is set largely by 
c, and depends only weakly upon temperature. How- 
ever, if c were strongly temperature-dependent, as has 
been suggested [l6j , we would observe an hysteresis that 
becomes more pronounced with increasing temperature. 

In Figure [5] we contrast the character of hysteresis ob- 
served in force-melting mode in the absence and presence 
of nicks. 
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FIG. 8: Illustration of the difference in the character of hys- 
teresis observed in force-melting mode when unpeeling is sup- 
pressed and allowed. Left panel: un-nicked A-DNA, c = 1.8, 
T = 25°C, 150 mM NaCl, pulling rates 1.8 and 3.6 /im/s. 
Melting driven hysteresis is asymmetric in character, and 
driven by the loop entropy exponent c. Equilibrium data 
(symbols) demonstrate that stretching is in equilibrium. For 
the slower pulling rate we show the percentage of molten 
DNA during stretching and shortening (top panel; arrows 
correspond to arrows on force-extension curve). Right set: 
once-nicked A-DNA, c = 1.8, T = 25°C, 150 mM NaCl, 
pulling rates 1.8, 3.6 and 14 jim/s. Unpeeling-driven hys- 
teresis is symmetric in character: stretching and shortening 
force-extension lines lie equally far from equilibrium (sym- 
bols). We observe a pulling rate-dependent overstretching 
force. For the slowest pulling rate we show the percentage of 
unpeeled DNA during stretching and shortening (top panel; 
arrows correspond to arrows on force-extension curve). 



B. B-to-S mode 



In B-to-S mode, when unpeeling is suppressed, we ob- 
serve a melting-driven hysteresis and variation in plateau 
height (Figure ■ These features are also observed in 
force-melting mode (Section IVI Ap . However, unlike in 
force-melting mode, for a fixed value of c this hysteresis 
is temperature dependent, rendered so by virtue of the 
temperature-dependent competition between S- and M- 
DNA. Melting-driven hysteresis sets in at higher temper- 
atures (or lower salt concentrations) than does unpeeling- 
generated hysteresis, because of the thermodynamic pre- 
dominance of U- over M-DNA. Note that the data shown 
in Figure [5] were obtained under conditions of 150 mM 
NaCl (a salt concentration frequently used in experi- 
ments) ; data in the presence of a nick, Figure 21 were 
obtained at a salt concentration of 500 mM (th e salt con- 
centration used in the experiments of Ref. [13|). 
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FIG. 9: B-to-S mode in the absence of nicks: Effect of temper- 
ature. Simulated force-extension data for un-nicked A-phage 
DNA at 150 mM NaCl and pulling rate 3000 nm/s. We show 
data for T = 21, 33, 37, 45, 55 and 65°C (offset left to right), 
with loop exponent c = 1.8, together with equilibrium data 
for temperatures 21°C and 45°C (symbols). At low temper- 
ature, overstretching is a transition from B- to S-DNA. In- 
creasing temperature causes melting, resulting in the onset 
of c-dependent hysteresis and a drop in plateau height. The 
overstretched state at high temperature is largely M-DNA. 



VII. CONCLUSIONS 

We have demonstrated that a simple model of the DNA 
overstretching phenomenon can furnish kinetic predic- 
tions on experimental length- and timescales, given as 
inputs simple experimental thermodynamic and elastic 
data. This approach offers the possibility of comparing in 
a detailed manner the two pictures of overstretching. One 
of our key results is the observation that hysteresis gen- 
erated by strand unpeeling has a character that depends 
crucially on the presence of the S-state of DNA, display- 
ing the asymmetry of experimentally-observed hysteresis 
only when S-DNA is allowed. We observe a progressive 
increase with temperature in the degree of this hystere- 
sis, owing to a competition between the entropy of U- 
DNA and the basepairing energy of S-DNA. We have 
also demonstrated that the nonextensive thermodynam- 
ics of ssDNA (loop entropy) is associated with an anoma- 
lous kinetics, which may be detectable in overstretching 
experiments performed upon un-nicked but torsionally 
unconstrained DNA. 

Our interpretation of force-melting theory yields kinet- 
ics in generally poor agreement with experiment. The 
chief cause of this disagreement is the instability of 
molten DNA to unpeeling: in force-melting mode, our 
model exhibits complete unpeeling at high forces, con- 
ferring upon the overstretching transition a symmetric 
hysteresis at all temperatures. This instability can be 
curbed within our model only by interposing between M- 
and U states a highly collective free energetic or kinetic 
barrier. 



However, we do not discount the possibility that a 
modified model of force-melting theory might better de- 
scribe experiment. Based on our work, we pose the fol- 
lowing questions. Is there a collective kinetic or thermo- 
dynamic effect that prevents unpeeling from nicks in the 
backbone, rendering melting stable or metastable relative 
to unpeeling? If so, might a high-force unpeeling follow 
a lower-force melting, thereby reproducing the two-stage 
transition shown in Figure [5]? Does tension- melted DNA 
differ fundamentally from thermally-melted DNA (other 
than in their respective lengths)? Does the loop factor of 
molten DNA change under tension? Can the loop factor 
be strongly suppressed without losing much of the ex- 
tensive entropy of melting? Our interpretation of force- 
melting theory assumes or results in simple answers to 
these questions: 'no' in each case. 

If some or all of these questions could be answered in 
the positive, we would expect better agreement between 
model predictions and experimental kinetics. Were melt- 
ing somehow stable or metastable to unpeeling (perhaps 
rendered so by hydrodynamic effects or effects of sec- 
ondary structure) we could in principle observe a high- 
force unpeeling transition following an initial melting, 
giving rise to force-extension data analagous to those pre- 
sented in Figure If the degree of this metastability 
was temperature-dependent, or if c were temperature- 
dependent, we might again capture the trend of increas- 
ingly prominent hysteresis with increasing temperature, 
such as shown in Figure [U Our results indicate that a 
change in c by roughly a factor of 2 over a 40° C range 
(i.e. c « 1 or 1.5 at 10°C and c « 2 at 50°C) should be 
sufficient to generate the degree of variation in hysteresis 
seen in experiment. Interestin gly, the variation of c with 
temperature suggested in Ref. [161 ] is consistent with this 
requirement. 

Determining the value of c under force would allow 
us to make detailed kinetic predictions when large-scale 
melting proliferates in strained DNA. To measure c, one 
could use the method proposed in Ref. (3f| . In this work 
the authors studied, theoretically, certain autocorrelation 
functions associated with basepair openings during DNA 
melting. They demonstrated that these quantities pos- 
sess a long-time behavior that is a simple function of 
c, and suggested that one could extract c by measuring 
these correlation functions using fluoresence correlation 
spectroscopy. We suggest that performing such measure- 
ments with the DNA in question held at different tensions 
may provide a means of discerning the force dependence 
of c. 

Measurement of the value of c under tension would 
allow our model to probe the kinetic behavior of large 
molten bubbles that develop in strained DNA. For a 
range of c (w 1 to 2) we expect there to be an anomalous 
kinetics associated with these bubbles. This may have 
significance for the behavior of the molecule in biological 
processes. 
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VIII. APPENDIX A 

Approximate entropy of loops under tension. The en- 
tropy associated with molten DNA has been shown to 
be important during overstretching from the strong tem- 
perature dependence of the overstretching force [f| ; from 
theoretical fits to experimental data (3lj ; and from molec- 
ular dynamics simulations [2(|. Associated with the ex- 
tensive entropy of melting is a nonextensive penalty for 
closing a loop, which to a first approximation may be 
calculated in the following manner [lj| . 

Consider an ideal random walk of n steps on a d- 
dimensional hypercubic lattice, constrained to return to 
its origin. Because the loop closes, the number of steps 
in the + and — directions in each dimension must be 
equal; the number of ways of ordering m such steps in 
one dimension is W(m) = m!/[(m/2)!] 2 . Consider a walk 
in which rn < n steps are taken in a favoured dimen- 
sion, corresponding to the direction of pulling. The case 
r = l/d is an unbiased walk. If the remaining (1 — r)n 
steps are distributed equally amongst the other d — 1 
dimensions, then the entropy of the walk is 

S(n, r) = k B In W(rn) + (d - l)k B In W{fn). (17) 

We have defined f = (1 — r)/(d — 1). The first term 
in Equation (|17p accounts for the number of ways of or- 
dering the steps in the 'pulling' dimension; the second 
accounts for the steps in the other d — 1 dimensions. 
Equation (|17|) can be approximated by Stirling's formula 
as 

S(n, r) w k B n In 2 — k B ca lnn r . (18) 

The first term in Equation (|18|) is the extensive entropy 
accounting for the 2 choices (+/— ) for each of the n 
steps of the walk. The second term is the nonextensive 
penalty for loop closure. The dependence upon aspect 
ratio r is confined to the effective loop length parameter 
n r = \irr 1 l d r 1 ~ 1 l d n. The exponent Co = d/2 does not 
depend on r. Under zero tension, in d — 3, this value is 
increased to ~ 1.8 by the constraint of self-avoidance [20], 
and further to ~ 2.1 if loop interactions are taken into 
account (2l| . 

Stretching indeed changes the nonextensive entropy of 
a loop. However, in the approximation that we con- 
sider here this change in entropy is independent of the 
length of the loop, and can be regarded as a constant 
(n-independent) correction for all loops. Such is also the 
case within the freely-jointed chain model (again ignor- 
ing excluded volume) in the constant force ensemble [23| . 



This calculation yields the entropic correction in terms 
of the pulling force /: 

ASioo P (/, r) « -k B c In (n) + k B In <?(/), (19) 

where 

^ = «/,)(«« + co t h(/ t >) . (20) 

Here / r = fP ss /(k B T) is a reduced force, P ss w 0.7 
nm is the persistence length of ssDNA, and C(x) = 
1 — [coth(x) — 1/x} 2 . The correction g(f) provides an 
approximate means of accounting for the change in loop 
entropy with force. Note that g(f) does not tend ex- 
actly to unity in the limit / — > 0, and so we regard this 
correction only as a rough estimate of the change in en- 
tropy of a loop under pulling. The emergent kinetics of 
our model does not depend on whether g{f) is present or 
absent; this term serves only to renormalize the surface 
tension of dsDNA-ssDNA interfaces (from ~ 11.5/cbT at 
zero tension to ~ 6fceT at 65 pN). However, when vol- 
ume exclusion is accounted for, the correction factor g{f) 
acquires a dependence upon loop length (23|. Instead of 
adopting this correction explicitly, we choose to vary c in 
order to probe the regime in which melting is appreciably 
hysteretic. 

IX. APPENDIX B 

What is the appropriate base-closing timescale for our 
coarse-grained model of DNA ?. As discussed in Section 
ITVl we use as the fundamental 'basepair closing' timescale 
of our model a value of 28 /j,s, derived from fluorescence 
correlation spectroscopy data [Tt], HH. In Ref. @ an 
NMR-motivated timescale of 10~ 8 s was used with a sim- 
ple model to derive kinetic predictions for the degree and 
nature of hysteresis of overstretching when unpeeling is 
permitted. Given that these kinetic predictions appear 
similar in nature to some of the data presented in this 
paper, we must explain why this can be so when two 
dramatically different timescales have been employed. 

We believe that the difference in timescales is not evi- 
dent because, chiefly, of the different dynamical protocols 
employed in the two papers. The dissipative dynamics of 
Ref Q is not equivalent to the Monte Carlo dynamics 
we use here. The former describes a dynamics that is 
dominated by the local free energy landscape of the in- 
homogeneous sequence, as we show below. While our 
Monte Carlo dynamics is indeed influenced by DNA se- 
quence, many of the features we describe in this paper, 
such as the asymmetry of hysteresis in the presence of 
S-DNA, do not depend solely on sequence heterogenity. 

In Figure [TO] we show data obtained from a model 
closely related to the kinetic model of refence 0. Our 
simple model consists of a 41 kbp fragment of A-DNA, 
with each basepair allowed to assume only two states. 
The first state is a hybrid 'B-S' phase, possessing an 
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elastic behavior that mirrors the form of the overstretch- 
ing plateau (see Ref. Q), and having the basepairing- 
stacking energy of the B-DNA of our Monte Carlo model. 
The second phase is the unpeeled phase, having the ther- 
modynamic and elastic behavior of the U state of our 
Monte Carlo model. A nick is placed between sites 2 and 
3. Starting from the fully hybridized state, we increment 
the length of the DNA at 1000 nm/s, according to a fun- 
damental timescale Tt, until the tension of the molecule 
reaches 140 pN. We then reverse the loading rate. We do 
not model the optical trap. 

Because only the 'B-S' and U states are accessible, the 
only relevant dynamical variable is the junction between 
these phases, which we call the unpeeled 'front'. We allow 
the front to advance or retreat along the inhomogeneous 
sequence according to a 'zero-temperature' dynamics: if 
a change in position of the front reduces the free energy 
of the system then this change is accepted; if not, it is re- 
jected. This dynamics is analogous to the dissipative dy- 
namics employed in Ref Q . We find that the kinetic pre- 
dictions of this model are dominated by sequence hetero- 
geneity. In Figure[TU] (top panel) we demonstrate that the 
nature (extent and 'roughness') of hysteresis does not de- 
pend strongly on the fundamental timescale Tf , provided 
that the nicked front has time to reach the local minimum 
of the free energy landscape. Whether this is so depends 



strongly on the nature of the landscape, which is influ- 
enced by any 'smoothing' that one might apply (in order 
to account on small lengthscales for thermally-driven be- 
havior Q). In Figure [TU1 (bottom panel) we show that 
averaging the local sequence using windows of 0, 10, 30 
and 100 bp influences strongly the zero-temperature dy- 
namics. 

The 'roughness' of unpeeling-generated hysteresis in 
our Monte Carlo model is not wholly determined by se- 
quence, but owes its form in part to the response of the 
optical trap and to thermal noise. Indeed, the 'sawteeth' 
features seen during relaxation differ from run to run; 
this is not the case for the dissipative dynamics described 
in this Appendix. Our conclusion from these data is 
that there is no direct contradiction between the basic 
timescale we have used in this paper (28 /is), in concert 
with Glauber dynamics, and the 10 -8 s employed as the 
integration step of the dissipative dynamics of Ref. H. 
These different fundamental timescales yield apparently 
similar predictions for the character of unpeeling-driven 
hysteresis because predictions are made using different 
dynamical protocols. Within our model, FCS timescales 
yield much better agreement with experiment when con- 
sidering the slow, high-force unpeeling transition follow- 
ing the initial overstretching plateau (Refs. [3, EH and 
FigureE]). 
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